Gas detection systems and methods with search directions

ABSTRACT

In some embodiments, at least one processor determines a range of search directions to a suspected gas leak source by receiving data representative of a plurality of wind direction measurements relative to a geo-referenced location of at least one gas concentration measurement point. The range of search directions extending from the location of the gas concentration measurement point is calculated according to a variability of the wind direction measurements.

CROSS REFERENCE TO RELATED APPLICATIONS

This application claims the benefit of U.S. provisional patentapplication 61/646,487 filed on May 14, 2012, titled “Gas DetectionSystems and Methods”, which application is hereby incorporated byreference in its entirety.

BACKGROUND

The invention relates to systems and methods for detecting gas leakssuch as a methane gas leaks.

A common means of distributing energy around the world is by thetransmission of gas, usually natural gas. In some areas of the worldmanufactured gasses are also transmitted for use in homes and factories.Gas is typically transmitted through underground pipelines havingbranches that extend into homes and other buildings for use in providingenergy for space and water heating. Many thousands of miles of gaspipeline exist in virtually every major populated area. Since gas ishighly combustible, gas leakage is a serious safety concern. Recently,there have been reports of serious fires or explosions caused by leakageof gas in the United States as the pipeline infrastructure becomesolder. For this reason, much effort has been made to provideinstrumentation for detecting small amounts of gas so that leaks can belocated to permit repairs.

Conventionally, search teams are equipped with gas detectors to locate agas leak in the immediate proximity of the detector. When the plume ofgas from a leak is detected, the engineers may walk to scan the areaslowly and in all directions by trial and error to find the source ofthe gas leak. This process may be further complicated by wind thatquickly disperses the gas plume. Such a search method is time consumingand often unreliable, because the engineer walks around with little orno guidance while trying to find the source of the gas leak.

Another approach to gas leak detection is to mount a gas leak detectioninstrument on a moving vehicle, e.g., as considered in U.S. Pat. No.5,946,095. A natural gas detector apparatus is mounted to the vehicle sothat the vehicle transports the detector apparatus over an area ofinterest at speeds of up to 20 miles per hour. The apparatus is arrangedsuch that natural gas intercepts a beam path and absorbs representativewavelengths of a light beam. A receiver section receives a portion ofthe light beam onto an electro-optical etalon for detecting the gas.Although a moving vehicle may cover more ground than a surveyor on foot,there is still the problem of locating the gas leak source (e.g., abroken pipe) if a plume of gas is detected from the vehicle. Thus, thereis still a need to provide a method and apparatus to locate the sourceof a gas leak quickly and reliably.

SUMMARY

According to one aspect, a method comprises employing at least oneprocessor to determine a range of search directions to a suspected gasleak source by receiving data representative of a plurality of winddirection measurements relative to a geo-referenced location of at leastone gas concentration measurement point. The range of search directionsextending from the location of the gas concentration measurement pointis calculated according to a variability of the wind directionmeasurements.

According to another aspect, an apparatus comprises at least oneprocessor programmed to determine a range of search directions to asuspected gas leak source by receiving data representative of aplurality of wind direction measurements relative to a geo-referencedlocation of at least one gas concentration measurement point. The rangeof search directions extending from the location of the gasconcentration measurement point is calculated according to a variabilityof the wind direction measurements.

According to another aspect, a non-transitory computer-readable mediumencodes instructions which, when executed by a computer system, causethe computer system to determine a range of search directions to asuspected gas leak source by receiving data representative of aplurality of wind direction measurements relative to a geo-referencedlocation of at least one gas concentration measurement point. The rangeof search directions extending from the location of the gasconcentration measurement point is calculated according to a variabilityof the wind direction measurements.

According to another aspect, a method comprises employing at least oneprocessor to determine a probable direction and an estimated maximumdistance to a suspected gas leak source by receiving data representativeof wind direction and wind speed relative to a geo-referenced locationof at least one gas concentration measurement point. The probabledirection to the suspected gas leak source relative to the location ofthe gas concentration measurement point is determined according to thedata representative of wind direction. The method also includes the stepof determining a maximum detection distance value representative of theestimated maximum distance from the suspected gas leak source at which aleak can be detected. The maximum detection distance value is determinedaccording to the data representative of wind speed.

According to another aspect, an apparatus comprises at least oneprocessor programmed to determine a probable direction and an estimatedmaximum distance to a suspected gas leak source by receiving datarepresentative of wind direction and wind speed relative to ageo-referenced location of at least one gas concentration measurementpoint. The probable direction to the suspected gas leak source relativeto the location of the gas concentration measurement point is determinedaccording to the data representative of wind direction. The processor isalso programmed to determine a maximum detection distance valuerepresentative of the estimated maximum distance from the suspected gasleak source at which a leak can be detected. The maximum detectiondistance value is determined according to the data representative ofwind speed.

According to another aspect, a non-transitory computer-readable mediumencodes instructions which, when executed by a computer system, causethe computer system to determine a probable direction and an estimatedmaximum distance to a suspected gas leak source by receiving datarepresentative of wind direction and wind speed relative to ageo-referenced location of at least one gas concentration measurementpoint. The probable direction to the suspected gas leak source relativeto the location of the gas concentration measurement point is determinedaccording to the data representative of wind direction. Also determinedis a maximum detection distance value representative of the estimatedmaximum distance from the suspected gas leak source at which a gas leakcan be detected. The maximum detection distance value is determinedaccording to the data representative of wind speed.

BRIEF DESCRIPTION OF THE DRAWINGS

The foregoing aspects and advantages of the present invention willbecome better understood upon reading the following detailed descriptionand upon reference to the drawings where:

FIG. 1 shows a gas leak detection apparatus according to someembodiments of the present invention.

FIG. 2 illustrates hardware components of a computer system according tosome embodiments of the present invention.

FIG. 3 shows a number of application or modules running on a clientcomputer system and a corresponding server computer system according tosome embodiments of the present invention.

FIG. 4 is a schematic drawing of a screen shot on a graphical userinterface displaying survey results on a street map according to someembodiments of the present invention.

FIG. 5 is a schematic drawing of a screen shot on a graphical userinterface with GPS indicators according to some embodiments of thepresent invention.

FIG. 6 is a schematic drawing of a screen shot on a graphical userinterface with weather station status indicators according to someembodiments of the present invention.

FIG. 7 is a schematic drawing of a screen shot on a graphical userinterface with map controls according to some embodiments of the presentinvention.

FIG. 8 is a schematic diagram of three search area indicators accordingto some embodiments of the present invention.

FIG. 9 is a schematic diagram illustrating wind lines relative to thepath of a mobile gas measurement device for detecting or not detecting agas leak from a potential gas leak source according to some embodimentsof the present invention.

FIG. 10 is a schematic diagram of wind direction and a path of a mobilegas measurement device used to estimate a probability of detection of agas leak from a potential gas leak source at one or more measurementpoints along the path according to some embodiments of the presentinvention.

FIG. 11 is a graph of probability density vs. wind directions forestimating a probability of detection of a gas leak from a potential gasleak source according to some embodiments of the present invention.

FIG. 12 is a flow chart showing steps for performing a gas leak surveyaccording to some embodiments of the present invention.

FIG. 13 is a flow chart showing steps for generating a search areaindicator according to some embodiments of the present invention.

FIG. 14 is a flow chart showing steps for calculating a boundary of asurvey area according to some embodiments of the present invention.

FIG. 15 is a flow chart showing steps for displaying layers overlaid orsuperimposed on a street map according to some embodiments of thepresent invention.

FIG. 16 is a graph of vertical dispersion coefficients of a gas plume asa function of downwind distance from a gas leak source according to someembodiments of the present invention.

FIG. 17 is a graph of crosswind dispersion coefficients of a gas plumeas a function of downwind distance from a gas leak source according tosome embodiments of the present invention.

FIG. 18 is a table of dispersion coefficients for various atmosphericconditions according to some embodiments of the present invention.

DETAILED DESCRIPTION

Apparatus and methods described herein may include or employ one or moreinterconnected computer systems such as servers, personal computersand/or mobile communication devices, each comprising one or moreprocessors and associated memory, storage, input and display devices.Such computer systems may run software implementing methods describedherein when executed on hardware. In the following description, it isunderstood that all recited connections between structures can be directoperative connections or indirect operative connections throughintermediary structures. A set of elements includes one or moreelements. Any recitation of an element is understood to refer to atleast one element. A plurality of elements includes at least twoelements. Unless otherwise required, any described method steps need notbe necessarily performed in a particular illustrated order. A firstelement (e.g. data) derived from a second element encompasses a firstelement equal to the second element, as well as a first elementgenerated by processing the second element and optionally other data.Making a determination or decision according to a parameter encompassesmaking the determination or decision according to the parameter andoptionally according to other data. Unless otherwise specified, anindicator of some quantity/data may be the quantity/data itself, or anindicator different from the quantity/data itself. Computer programsdescribed in some embodiments of the present invention may bestand-alone software entities or sub-entities (e.g., subroutines, codeobjects) of other computer programs. Computer readable media encompassstorage (non-transitory) media such as magnetic, optic, andsemiconductor media (e.g. hard drives, optical disks, flash memory,DRAM), as well as communications links such as conductive cables andfiber optic links. According to some embodiments, the present inventionprovides, inter alia, computer systems programmed to perform the methodsdescribed herein, as well as computer-readable media encodinginstructions to perform the methods described herein.

FIG. 1 shows a gas leak detection system 10 according to someembodiments of the present invention. System 10 comprises a serviceprovider server computer system 12 and a set of client computer systems18 a-c all connected through a wide area network 16 such as theInternet. Client computer systems 18 a-c may be personal computers,laptops, smartphones, tablet computers and the like. A vehicle 24 suchas an automobile may be used to carry at least some client computersystems (e.g. an exemplary client computer system 18 c) and associatedhardware including a mobile gas measurement device 26, a location/GPSmeasurement device 30, and a wind measurement device 32. In a preferredembodiment, the mobile gas measurement device 26 may be a Picarroanalyzer using Wavelength-Scanned Cavity Ring Down Spectroscopy (CRDS),available from Picarro, Inc., Santa Clara, Calif. Such analyzers may becapable of detecting trace amounts of gases such as methane, acetylene,carbon monoxide, carbon dioxide, hydrogen sulfide, and/or water. Inparticular applications suited for detection of natural gas leaks, aPicarro G2203 analyzer capable of detecting methane concentrationvariations of 3 ppb may be used. Wind measurement device 32 may includea wind anemometer and a wind direction detector (e.g. wind vane). GPSmeasurement device 30 may be a stand-alone device or a device built intoclient computer system 18 c.

FIG. 2 schematically illustrates a plurality of hardware components thateach computer system 18 may include. Such computer systems may bedevices capable of web browsing and have access to remotely-hostedprotected websites, such as desktop, laptop, tablet computer devices, ormobile phones such as smartphones. In some embodiments, computer system18 comprises one or more processors 34, a memory unit 36, a set of inputdevices 40, a set of output devices 44, a set of storage devices 42, anda communication interface controller 46, all connected by a set of buses38. In some embodiments, processor 34 comprises a physical device, suchas a multi-core integrated circuit, configured to execute computationaland/or logical operations with a set of signals and/or data. In someembodiments, such logical operations are delivered to processor 34 inthe form of a sequence of processor instructions (e.g. machine code orother type of software). Memory unit 36 may comprise random-accessmemory (RAM) storing instructions and operands accessed and/or generatedby processor 34. Input devices 40 may include touch-sensitiveinterfaces, computer keyboards and mice, among others, allowing a userto introduce data and/or instructions into system 18. Output devices 44may include display devices such as monitors. In some embodiments, inputdevices 40 and output devices 44 may share a common piece of hardware,as in the case of touch-screen devices. Storage devices 42 includecomputer-readable media enabling the storage, reading, and writing ofsoftware instructions and/or data. Exemplary storage devices 42 includemagnetic and optical disks and flash memory devices, as well asremovable media such as CD and/or DVD disks and drives. Communicationinterface controller 46 enables system 18 to connect to a computernetwork and/or to other machines/computer systems. Typical communicationinterface controllers 46 include network adapters. Buses 38 collectivelyrepresent the plurality of system, peripheral, and chipset buses, and/orall other circuitry enabling the inter-communication of devices 34-46 ofcomputer system 18.

FIG. 3 shows a number of applications or modules running on an exemplaryclient computer system 18 and corresponding server computer system 12.Authentication applications 54, 62 are used to establish securecommunications between computer systems 12, 18, allowing client computersystem 18 selective access to the data of a particular customer or useraccount. A client data collection module 56 collects real-time gasconcentration, location data such as global positioning system (GPS)data, as well as wind speed and wind direction data. A graphical userinterface (GUI) module 58 is used to receive user input and displaysurvey results and other GUI displays to system users. A client-sidereal-time data processing module 60 may be used to perform at least someof the data processing described herein to generate survey results frominput data. Data processing may also be performed by a server-side dataprocessing module 68. Server computer system 12 also maintains one ormore application modules and/or associated data structures storingcurrent and past survey results 64, as well as application modulesand/or data structures storing reference data 66 such as platsindicating the geographic locations of natural gas pipelines.

FIG. 4 is a schematic drawing of a screen shot on a graphical userinterface, displaying survey results on a street map 70 according tosome embodiments of the present invention. The GUI screenshots are mostpreferably displayed on a client device in the vehicle, which may beconnected to a server as described above. The illustrated screenshotsshow both exemplary user input, which may be used to control systemoperation, and exemplary real-time displays of collected/processed data.In the example, it includes the geo-referenced street map 70 showingplat lines 72. The plat lines 72 are preferably derived from gas companyrecords. An active pipeline plat boundary 71 may also be displayed onthe map 70. A user-selectable button 96 may be selected to overlay aselected pipeline plat on the map 70. Superimposed on the map 70 are oneor more lines (preferably in a distinguishing color not shown in patentdrawings) indicating the path 74 driven by the vehicle with the mobilegas measurement device on one or more gas survey routes. In thisexample, the path 74 shows the vehicle U-turned at the Y-shapedintersection. Optionally, a current location icon 75 may be overlaid onthe map 70 to indicate the current surveyor location, e.g., the positionof the vehicle with a gas measurement device and wind measurementdevice. A user-selectable button 94 may be selected to center the map 70by current surveyor location. Also provided is a user-selectable startbutton 102 and stop button 100 to start/stop capturing gas for analysis.An analyzer control button 104 is user-selectable to control analyzeroperations (e.g., shut down, start new trace, run isotopic analysis,etc.).

Peak markers 76 show the locations along the path 74 where peaks in thegas concentration measurements, which satisfy the conditions for beinglikely gas leak indications, were identified. The colors of the peakmarkers 76 may be used to distinguish data collected on different runs.The annotations within the peak markers 76 show the peak concentrationof methane at the locations of those measurement points (e.g., 3.0, 2.6,and 2.0 parts per million). An isotopic ratio marker 77 may be overlaidon the map 70 to indicate isotopic ratio analysis output and tolerance(e.g., −34.3+/−2.2). Also displayed on the map 70 are search areaindicators 78, preferably shown as a sector of a circle having adistinguishing color. Each of the search area indicators 78 indicates asearch area suspected to have a gas leak source. The opening angle ofthe search area indicator 78 depicts the variability in the winddirection. The axis of the search area indicator 78 (preferably an axisof symmetry) indicates the likely direction to the potential gas leaksource. Also displayed on the map 70 are one or more survey areaindicators 80 (shown as hatched regions in FIG. 4) that indicate asurvey area for a potential gas leak source. The survey area indicator80 adjoins the path 74 and extends in a substantially upwind directionfrom the path. The survey area marked by each indicator 80 is preferablydisplayed as a colored swath overlaid or superimposed on the map 70. Forexample, the colored swaths may be displayed in orange and green for tworuns. In preferred embodiments, the parameters of the search areaindicators 78 and the survey area indicators 80 (described in greaterdetail with reference to FIGS. 8-11 below) are derived from measurementsof the wind, the velocity of the vehicle, and optionally the prevailingatmospheric stability conditions.

Referring still to FIG. 4, the surveyor user interface also preferablyincludes a real-time CH4 concentration reading 82. A wind indicatorsymbol 84 preferably displays real-time wind information, which may becorrected for the velocity vector of the vehicle to represent the truewind rather than the apparent wind when the vehicle is in motion.Average wind direction is preferably indicated by the direction of thearrow of the wind indicator symbol 84, while wind direction variabilityis indicated by the degree of open angle of the wedge extending from thebottom of the arrow. Wind speed is preferably indicated by the length ofthe arrow in the wind indicator symbol 84. An internet connectionindicator 98 blinks when the internet connection is good. A datatransfer status button 86 is user-selectable to display data transferstatus (e.g., data transfer successful, intermittent data transfer, ordata transfer failed). An analyzer status button 88 is user-selectableto display current analyzer status such as cavity pressure, cavitytemperature, and warm box temperature. A map control button 106 isuser-selectable to open a map controls window with user-selectable layeroptions, discussed below with reference to FIG. 7.

FIG. 5 is a schematic drawing of a screen shot on the graphical userinterface, displaying a GPS status window 91, according to someembodiments of the present invention. A user-selectable GPS statusbutton 90 may be selected to open the GPS status window 91. The GPSstatus window 91 preferably includes indicators of the current GPSstatus, such as “GPS OK”, “Unreliable GPS signal”, “GPS Failed”, or “GPSNot Detectable”.

FIG. 6 is a schematic drawing of a screen shot on the graphical userinterface, displaying a weather station status window 93, according tosome embodiments of the present invention. A user-selectable weatherstation status button 92 may be selected to open the weather stationstatus window 93. The weather station status window 93 preferablyincludes indicators of the current weather station status, such as“Weather Station OK”, “Weather Station Failed”, or “Weather Station NotDetectable”. Weather station data are preferably received in real-timeand may include wind data and atmospheric stability conditions datarelevant to the area being surveyed.

FIG. 7 is a schematic drawing of a screen shot on the graphical userinterface, displaying a map control window 95. Various elementsdisplayed on the map 70 are regarded as layers which may be turned on oroff. In this example, map controls window 95 includes sixuser-selectable buttons named “Hide Peak Markers”, “Hide Search AreaIndicators”, “Minimum Amplitude”, “Hide Isotopic Analysis”, “Hide Fieldof View”, and “Hide Plat Outline”. The “Hide Peak Markers” button may beselected so that the markers indicating peak gas concentrationmeasurements are not displayed on the map 70. The “Hide Search AreaIndicators” button may be selected so that the search area indicatorsare not displayed on the map 70. The “Minimum Amplitude” button may beselected so that gas concentration peaks not meeting a minimum amplituderequirement are not displayed on the map 70. The “Hide IsotopicAnalysis” button may be selected so that isotopic ratio analysisinformation is not displayed on the map 70 next to the peak markers. The“Hide Field of View” button may be selected so that the survey areaindicator(s) are not displayed on the map 70. The “Hide Plat Outline”button may be selected so that the plat lines are not displayed on themap 70.

FIG. 8 is a schematic diagram of three search area indicators 78 a, 78b, and 78 c according to some embodiments of the present invention. Eachof the search area indicators 78 a, 78 b, and 78 c has a respective axis108 a, 108 b, and 108 c indicating a representative wind directionrelative to a geo-referenced location of a corresponding gasconcentration measurement point M1, M2, and M3. The gas concentrationmeasurement points M1, M2, and M3 are positioned along the path 74traveled by the vehicle 24 that carries a GPS device, a mobile gasmeasurement device, and wind measurement device for taking winddirection measurements and wind speed measurements. Each of the searcharea indicators, such as the search area indicator 78 c, preferably hasa width W relative to its axis 108 c. The width W is indicative of awind direction variability associated with wind direction measurementsin the area of the gas concentration measurement point M3. In preferredembodiments, the width W is indicative of a variance or standarddeviation of the wind direction measurements. Also in preferredembodiments, the search area indicator 78 c has the shape of a sector ofa circle, with the center of the circle positioned on the map at thelocation of the gas concentration measurement point M3. Most preferably,the angle A subtended by the sector of the circle is proportional to astandard deviation of the wind direction measurements taken at or nearbythe measurement point M3. For example, the angle A may be set to a valuethat is twice the angular standard deviation of the wind directionmeasurements. It is not necessary to display the gas concentrationmeasurement points M1, M2, and M3 on the map along with the search areaindicators 78 a, 78 b, and 78 c. As previously shown in FIGS. 4 and 7,the measurement points and associated gas concentration measurements arepreferably map layer options for an end-user that may be turned on oroff.

Referring again to FIG. 8, the axis 108 c of the search area indicator78 c is preferably an axis of symmetry and points in a representativewind direction relative to the gas concentration measurement point M3.The representative wind direction is preferably a mean, median or modeof the wind direction measurements taken at or nearby the measurementpoint M3, and indicates the likely direction to a potential gas leaksource. The wind direction measurements may be taken from the vehicle 24as it moves and converted to wind direction values relative to theground (e.g., by subtracting or correcting for the velocity vector ofthe vehicle). In some embodiments, the axis 108 c has a length Lindicative of a maximum detection distance value representative of anestimated maximum distance from a potential gas leak source at which agas leak from the source can be detected. For example, the length may beproportional to the maximum detection distance value, or proportional toa monotonically increasing function of the maximum detection distancevalue, such that longer maximum detection distance values arerepresented by longer axis lengths. In preferred embodiments, themaximum detection distance value and corresponding length L aredetermined according to data representative of wind speed in the searcharea. In some embodiments, the maximum detection distance value and thecorresponding length L are determined according to data representativeof atmospheric stability conditions in the search area. Each of thesearch area indicators 78 a, 78 b, and 78 c may thus provide a visualindication of a likely direction and estimated distance to a potentialgas leak source. Although a sector of a circle is the presentlypreferred shape for a search area indicator, alternative shapes for asearch area indicator include, but are not limited to, a triangle, atrapezoid, or a wedge.

FIG. 9 is a schematic diagram illustrating an example of detecting ornot detecting a gas leak from a potential gas leak source, according tosome embodiments of the present invention. An indicator of a survey area(also sometimes referred to as a “field of view”) is intended as anindication of how well the measurement process surveys the area aroundthe path 74 traveled by the vehicle 24 that carries a GPS device, amobile gas measurement device, and wind measurement device. The surveyarea indicator is designed such that if a potential gas leak source islocated in the survey area and has a rate of leakage meeting a minimumleak rate condition, then an estimated probability of detection of a gasleak from the potential gas leak source at one or more measurementpoints P along the path 74 satisfies a probability condition.

Whether or not a potential gas leak source of a given strength isdetectable by a gas measurement device of a given sensitivity depends onthe separation distance of the source from the gas measurement deviceand on whether the wind is sufficient to transport gas from the gas leaksource to the gas measurement device at some point along the path 74. Insome embodiments, a physical model is employed that relates the measuredgas concentration peak at the location of the vehicle 24 (in ppm, forexample) to the emission rate of the potential gas leak source (ing/sec, for example) and the distance between the source and thedetection point.

There are multiple possible models that describe the propagation of agas leak as a plume through the atmosphere. One well-validated physicalmodel for a plume (Gifford, F. A., 1959. “Statistical properties of afluctuating plume dispersion model”. Adv. Geophys, 6, 117-137) is tomodel the plume as a Gaussian distribution in the spatial dimensionstransverse to the wind direction, or (for a ground level source), theconcentration c (x, y, z) at a distance x downwind, y crosswind, and ata height z from a gas leak source of strength Q located on the ground isgiven by Equation (1):

$\begin{matrix}{{C\left( {x,y,z} \right)} = {\frac{Q}{\pi\; v\;\sigma_{y}\sigma_{z}}e^{{{{- y^{2}}/2}\;\sigma_{y}^{2}} - {{z^{2}/2}\sigma_{z}^{2}}}}} & \lbrack 1\rbrack\end{matrix}$where ν is the speed of the wind, and the plume dispersion half-widthsσ_(y) and σ_(z) depend on x via functions that are empiricallydetermined for various atmospheric stability conditions.

If we consider the plume center, where y=z=0, the concentration at thecenter is given by Equation (2):

$\begin{matrix}{C_{peak} = \frac{Q}{\pi\; v\;\sigma_{y}\sigma_{z}}} & \lbrack 2\rbrack\end{matrix}$

The dimensions of the Gaussian distribution horizontally and vertically,half-widths σ_(y) and σ_(z), increase with increasing distance from thesource. The amount they increase can be estimated from measurements ofwind speed, solar irradiation, ground albedo, humidity, and terrain andobstacles, all of which influence the turbulent mixing of theatmosphere. However, if one is willing to tolerate somewhat moreuncertainty in the distance estimation, the turbulent mixing of theatmosphere can be estimated simply from the wind speed, the time of day,and the degree of cloudiness, all of which are parameters that areavailable either on the vehicle 24 or from public weather databases inreal time. Using these available data, estimates of the Gaussian widthparameters can be estimated using the Pasquill-Gifford-Turner turbulencetyping scheme (Turner, D. B. (1970). “Workbook of atmospheric dispersionestimates”. US Department of Health, Education, and Welfare, NationalCenter for Air Pollution Control), or modified versions of this scheme.

For a given sensitivity of the gas measurement device, there is aminimum concentration which may be detected. Given a gas leak source ofstrength greater than or equal to the minimum concentration, the sourcewill be detected if it is closer than an estimated maximum distanceX_(max), where this is the distance such that σ_(y)σ_(z)=Q/(π νc). Ifthe wind is blowing gas directly from the gas leak source to the gasmeasurement device, the estimated maximum distance X_(max) is thedistance beyond which the source may be missed. This estimated maximumdetection distance may depend upon atmospheric stability conditions aswell as wind speed. The formula diverges to infinity when the wind speedis very small, so it is advisable to set a lower limit (e.g., 0.5 m/s)for this quantity.

The minimum leak rate Q_(min) is determined by the requirements of theapplication. For natural gas distribution systems, a minimum leak rateof 0.5 scfh (standard cubic feet per hour) may be used; below thislevel, the leak may be considered unimportant. Other minimum leaks rates(e.g. 0.1 scfh, 1 scfh, or other values within or outside this range)may be used for natural gas or other leak detection applications. Theminimum detection limit of the plume C_(min) is given either by the gasdetection instrument technology itself, or by the spatial variability ofmethane in the atmosphere when leaks are not present. A typical valuefor C_(min) is 30 ppb (parts-per-billion) above the background level(typically 1,800 ppb). Given these two values for Q_(min) and C_(min),and by predicting σ_(y) and σ_(z) given atmospheric measurements (orwith specific assumptions about the state of the atmosphere, such as thestability class), one may then determine the estimated maximum detectiondistance X_(max) by determining the value for X_(max) that satisfies thefollowing equality, Equation (3):

$\begin{matrix}{C_{m\; i\; n} = {\frac{Q_{m\; i\; n}}{\pi\; v\;\sigma_{y}\sigma_{z}}.}} & \lbrack 3\rbrack\end{matrix}$

In some embodiments the relationship between σ_(y) and σ_(z) and X_(max)is provided by a functional relationship, a lookup table, or similarmethod. Because σ_(y) and σ_(z) are monotonically increasing functionsof X_(max), a unique value can be determined from this process. Forexample, one useful functional form is a simple power law, where thecoefficients a, b, c, and d depend on atmospheric conditions:σ_(y)=ax^(b); σ_(z)=cx^(d).

In some embodiments, the concentration C measured close to the ground ofa Gaussian plume due to a gas leak source on the ground depends on therate of emission Q of the source, the distance x between the source andthe gas measurement device, and the speed of the wind blowing from thesource to the gas measurement device, in accordance with an expressionof the form (Equation 4):

$\begin{matrix}{C = \frac{Q}{\pi\; v\;{\sigma_{y}(x)}{\sigma_{z}(x)}}} & \lbrack 4\rbrack\end{matrix}$The expressions for σ_(y)(x) and σ_(z)(x) depend on the stability classof the atmosphere at the time of measurement. In some embodiments, thestability class of the atmosphere is inferred from the answers to a setof questions given to the operator, or from instruments of the vehicle,or from data received from public weather databases. As shown in thetable of FIG. 18, coefficients A, B, C, D, E and F may depend on surfacewind speed and atmospheric conditions such as day or night, incomingsolar radiation, and cloud cover. Mathematical forms for σ_(y)(x) andσ_(z)(x) are documented in Section 1.1.5 of the User's Guide forIndustrial Source Complex (ISC3), Dispersion Models Vol. 2 (USEnvironmental Protection Agency document EPA-454/B955-003b September1995). Given the sensitivity of the gas measurement device and the rateof emission of the smallest potential gas leak source of interest,equation (4) may be solved to find the estimated maximum distanceX_(max) beyond which a potential gas leak source may be missed by thegas measurement device.

FIG. 16 is a graph of vertical σ_(z)(x) dispersion coefficients of a gasplume as a function of downwind distance from a gas leak sourceaccording to some embodiments of the present invention. FIG. 17 is agraph of crosswind σ_(y)(x) dispersion coefficients of a gas plume as afunction of downwind distance from a gas leak source according to someembodiments of the present invention. The graphs are from de Nevers,2000, Air Pollution Control Engineering, The McGraw-Hill Companies, Inc.The dispersion coefficients are functions of downwind distance x. Inthis example, dispersion coefficients are calculated based onatmospheric stability. The table of FIG. 18 gives the atmosphericstability class as a function of wind speed, day or night, cloud cover,and solar radiation. In some embodiments, the dispersion coefficientsand/or the estimated maximum distance X_(max) may depend upon an urbanor rural environment for the gas concentration measurements and plumedispersion. For example, the estimated maximum distance X_(max) may beless in an urban environment with buildings or other structures than ina rural environment.

The actual distance at which a gas leak source may be detected isreduced if there is some variability or uncertainty in the direction ofthe wind. This is because there is a probability that the wind blows gasin a direction such that it does not intercept the path 74 of thevehicle 24 (FIG. 9). In practice this uncertainty is usually larger thanthe intrinsic angular uncertainty σ_(y)/x implied by the Gaussian plumemodel. In order to determine the effective survey area of the mobile gasmeasurement device, assume for this example that the wind speed remainsapproximately constant within a time interval −T<t<T bounding the timet=0 at which the vehicle 24 passes through a particular point P on thepath 74, but that the wind direction (angle) is distributed as aGaussian with a known mean and standard deviation.

As shown in FIG. 9, we consider the line 110 through the measurementpoint P pointing toward the direction of the mean wind, and whether acandidate point Q on this line qualifies to be within the boundary ofthe survey area (i.e., within the field of view of the mobile gasmeasurement device of the vehicle 24). We also consider drawing a samplefrom the distribution of wind directions and drawing a line through thecandidate point Q in this direction. If this line intersects the path 74of the vehicle 24 within the time interval −T<t<T, and the distance fromthe candidate point Q to the point of intersection with the path 74 isless than or equal to the estimated maximum distance X_(max), then thisis regarded as detectable by the mobile gas measurement device since thepotential gas leak source at the candidate point Q would have beendetected along the path 74. The quantity T sets the time interval duringwhich it is expected to detect the gas coming from the candidate point Qat measurement point P. Theoretically, the time interval can be large,but it may not be reasonable to assume that the wind statistics remainunchanged for an extended period of time. In some embodiments, the winddirection measurements are taken during a time interval less than orequal to about 2 minutes, during which time interval a gas concentrationis measured at the gas concentration measurement point P. Morepreferably, the time interval is in the range of 10 to 20 seconds.

FIG. 10 is a schematic diagram showing the estimation of a probabilityof detection at the measurement point P of a gas leak from a potentialgas leak source at the candidate point Q, according to some embodimentsof the present invention. The probability of detection at measurementpoint P is estimated according to an angle θ subtended by a segment 79of the path 74 relative to the candidate point Q for the potential gasleak source. The path segment 79 is positioned within a distance of thecandidate point Q that is less than or equal to the estimated maximumdistance X_(max). The probability of detection is preferably estimatedaccording to a cumulative probability of wind directions with respect tothe subtended angle θ. The cumulative probability of wind directions maybe determined according to a representative wind direction (e.g., amean, median, or mode of the wind direction measurements) and a winddirection variability (e.g., variance or standard deviation) calculatedfrom the wind direction measurements.

The candidate point Q is deemed to be within the boundary of the surveyarea if the probability of successful detection of a potential gas leaksource at the candidate point Q, over the distribution of winddirections, satisfies a probability condition. In some embodiments, theprobability condition to be satisfied is an estimated probability ofsuccessful detection greater than or equal to a threshold value,typically set at 70%. In general, as the candidate point Q is moved afarther distance from the gas concentration measurement point P, therange of successful angles becomes smaller and the probability ofsuccess decreases, reaching a probability threshold at the boundary ofthe territory deemed to be within the survey area.

FIG. 11 is a graph of probability density vs. wind directions forestimating a probability of detection of a gas leak from a potential gasleak source, according to some embodiments of the present invention. Thearea under the curve spans a range of possible angles θ for thesuccessful detection of a potential gas leak from a candidate point. Theprobability density is preferably generated as a Gaussian or similardistribution from the calculated mean and standard deviation of the winddirection measurements in the area of the gas concentration measurementpoint P, FIG. 10. If the angle θ subtended by the path segment 79relative to the candidate point Q encompasses a percentage of possiblewind vectors that is greater than equal to a threshold percentage (e.g.,70%, although the percentage may be adjusted to other values such as50%, 60%, 67%, 75%, 80%, or 90% in some embodiments), and if thedistance from the candidate point Q to the measurement point P is lessthan the estimated maximum distance X_(max), then the candidate point Qis deemed to be within the survey area.

The above process is repeated as different measurement points along thepath 74 are chosen and different candidate points are evaluated for theprobability of successful detection of a potential gas leak source. Thecumulative distribution of the wind direction function together with aroot finding algorithm are useful for efficiently determining theboundary of the survey area. For example, referring again to FIG. 10,the root finding algorithm may consider candidate points along the lineof mean wind direction starting at the estimated maximum distanceX_(max) from measurement point P, and iteratively (e.g. using abisection or other method) moving closer to the measurement point Palong the mean wind direction line until the angle θ subtended by thepath segment 79 is sufficient to meet the probability threshold, asdetermined from the cumulative probability of wind directions over thesubtended angle θ, FIG. 11. Referring again to FIG. 4, the survey areaindicator 80 may be displayed on the map 70 as a colored “swath”adjoining the path 74 and extending in a substantially upwind directionfrom the path.

FIG. 12 is a flow chart showing a sequence of steps to perform a gasleak survey according to some embodiments of the present invention. Instep 200, the survey program is started, preferably by an operator inthe vehicle using a graphical user interface (GUI). The operator beginsto drive the vehicle on a survey route while the GUI displays a streetmap (FIG. 4). Gas concentration measurements are preferably performedrapidly along the survey route (e.g., at a rate of 0.2 Hz or greater,more preferably 1 Hz or greater). This enables the practice of drivingthe vehicle at normal surface street speeds (e.g., 35 miles per hour)while accumulating useful data. The gas concentration is measuredinitially as a function of time, and is combined with the output of theGPS receiver in order to obtain the gas concentration as a function ofdistance or location. Interpolation can be used to sample the data on aregularly spaced collection of measurement points. The concentration ofmethane typically varies smoothly with position, for the most part beingequal to the worldwide background level of 1.8 parts per milliontogether with enhancements from large and relatively distant sourcessuch as landfills and marshes.

In step 210, at least one processor (e.g. of a client device, serverdevice, or a combination) receives data representative of measured gasconcentrations, wind direction measurements, wind speed measurements,and GPS data. In decision block 220, it is determined if a peak in gasconcentration is identified. A peak may be identified from a gasconcentration measurement above a certain threshold (or within a certainrange), or exceeding background levels by a certain amount, which may bepredetermined or user-selected. In some embodiments, the gasconcentration and GPS data are analyzed using a peak-location method,and then each identified peak is subsequently fit (using linear ornonlinear optimization) for center and width. The functional form usedfor this fitting step may be a Gaussian pulse, since a Gaussian iscommonly the expected functional form taken by gas plumes propagatingthrough the atmosphere.

If a peak in gas concentration is not identified, then the programproceeds to step 250. If a peak in gas concentration is identified, thena peak marker is generated in step 230. The peak marker may be displayedon the map as a user-selectable layer, as previously discussed withreference to FIG. 4. In step 240, a search area indicator is generatedto indicate the likely location of a gas leak source corresponding tothe identified peak in gas concentration. The search area indicator maybe displayed on the map as a user-selectable layer, as shown in FIG. 4.In step 250, the survey area boundary is calculated, and a survey areaindicator may be displayed on the map as a user-selectable layer(hatched region in FIG. 4). In decision step 260, it is determined ifthe operator wishes to continue surveying (e.g., by determining if the“Stop Survey” button has been selected). If yes, the survey programreturns to step 210. If not, the survey results are stored in memory instep 270 (e.g., in the survey results 64 of FIG. 3), and the surveyprogram ends.

FIG. 13 is a flow chart showing a sequence of steps performed togenerate a search area indicator according to some embodiments of thepresent invention. When a local enhancement in the gas concentration isdetected, the likely direction and estimated distance to the potentialgas leak source is preferably calculated from data representative ofwind direction and wind speed measured during a time interval just priorto or during which the gas concentration was measured. The time intervalis preferably fewer than 2 minutes, and more preferably in the range of5 to 20 seconds. Calculating statistics from wind measurements mayrequire some conversion if the measurements are made using sensors on amoving vehicle. A sonic anemometer is preferably used to measure windalong two perpendicular axes. Once the anemometer has been mounted tothe vehicle, these axes are preferably fixed with respect to thevehicle. In step 241, wind speed and wind direction values that weremeasured relative to the vehicle are converted to wind speed and winddirection values relative to the ground by subtracting the velocityvector of the vehicle, as obtained from the GPS data. When the vehicleis stationary, GPS velocity may be ineffective for determining theorientation of the vehicle and wind direction, so it is preferable touse a compass (calibrated for true north vs. magnetic north) in additionto the anemometer.

In step 242, wind statistics are calculated from the converted windvalues to provide the parameters for the search area indicator. Thestatistics include a representative wind direction that is preferably amean, median, or mode of the wind direction measurements. The statisticsalso include a wind direction variability, such as a standard deviationor variance of the wind direction measurements. In step 243, an angularrange of search directions, extending from the location of the gasconcentration measurement point where the local enhancement wasdetected, is calculated according to the variability of the winddirection measurements. In optional step 244, atmospheric conditionsdata are received. Step 245 is determining a maximum detection distancevalue representative of the estimated maximum distance from thesuspected gas leak source at which a leak can be detected. In someembodiments, the maximum detection distance value is determinedaccording to Equation (3) or Equation (4), and the data representativeof wind speed and/or atmospheric stability conditions. Alternatively,the maximum detection distance value may be a predetermined number, auser-defined value, empirically determined from experiments, or a valueobtained from a look-up table. In step 246, the search area indicator isgenerated with the determined parameters, previously discussed withreference to FIG. 8.

FIG. 14 is a flow chart showing a sequence of step performed tocalculate a boundary of a survey area according to some embodiments ofthe present invention. In step 251, wind speed and wind direction valuesthat were measured relative to the vehicle are converted to wind speedand wind direction values relative to the ground by subtracting thevelocity vector of the vehicle, as previously described in step 241above. In optional step 252, atmospheric conditions data are received.Step 253 is determining a maximum detection distance valuerepresentative of the estimated maximum distance from a suspected gasleak source at which a leak can be detected. In some embodiments, themaximum detection distance value is determined according to Equation (3)or Equation (4), and the data representative of wind speed and/oratmospheric stability conditions. Alternatively, the maximum detectiondistance value may be a predetermined number, a user-defined value,empirically determined from experiments, or a value obtained from alook-up table. In step 254, it is determined what angle θ is subtendedby a segment of the path of the vehicle relative to the candidate pointQ for the potential gas leak source. The path segment is positionedwithin a distance of the candidate point Q that is less than or equal tothe estimated maximum distance.

In step 255, a representative wind direction (e.g., a mean, median, ormode of the wind direction measurements) and a wind directionvariability (e.g., variance or standard deviation) are calculated fromthe wind direction measurements. In step 256, the probability ofdetection is estimated according to a cumulative probability of winddirections with respect to the subtended angle θ. In step 257, thesurvey area boundary is calculated with a probability threshold. Forexample, if the angle θ subtended by the path segment relative to thecandidate point encompasses a percentage of possible wind vectors thatis greater than equal to a threshold percentage (e.g., 70%), and if thedistance from the candidate point Q to the measurement point P is lessthan the estimated maximum distance X_(max), then the candidate point Qis deemed to be within the survey area. In decision step 258, it isdetermined if the survey area boundary function is to continue with thenext measurement point. If yes, steps 251-257 are repeated as differentmeasurement points along the path are chosen and different candidatepoints are evaluated for the probability of successful detection of apotential gas leak source. If not, then the boundary function ends.

FIG. 15 is a flow chart showing steps for displaying layers overlaid orsuperimposed on a street map according to some embodiments of thepresent invention. In step 310, a street map is displayed, preferably ona GUI visible to the operator in the vehicle. In step 320, the path ofthe vehicle with the mobile gas measurement device is displayed on themap. Various elements displayed on the map are regarded as layers whichmay be turned on or off. In this example, the map controls window (FIG.7) includes six user-selectable buttons named “Hide Peak Markers”, “HideSearch Area Indicators”, “Minimum Amplitude”, “Hide Isotopic Analysis”,“Hide Field of View”, and “Hide Plat Outline”. In decision step 330, itis determined if one or more of these layers is selected. If yes, theselected layer is displayed overlaid or superimposed on the street mapin step 340. If not, it is determined if the survey is to continue. Ifyes, display steps 310-350 are repeated. If not, the display options mayend.

The exemplary systems and methods described above allow a surveyor tolocate potential gas leak sources efficiently and effectively in highlypopulated areas. The search area indicators provide likely direction andestimated maximum distance to the source of detected gas leaks, whilethe survey area indicators provide an estimated statistical measure ofconfidence that an area was successfully surveyed for potential gasleaks. These aspects provide significant improvement in finding gas leaksources over conventional methods where engineers scan the area veryslowly and in all directions by trial and error to find the source of agas leak. These aspects also account for wind that may quickly dispersea gas plume.

It will be clear to one skilled in the art that the above embodimentsmay be altered in many ways without departing from the scope of theinvention. For example, gas leaks may include, but are not limited to:leaks from gas pipes or transportation systems (e.g., natural gasleaks), leaks from gas processing or handling facilities, and emissionsfrom gas sources into the environment (e.g., pollution, gas emissionfrom landfills, etc.). Gas concentration measurements are preferablyperformed rapidly (e.g., at a rate of 0.2 Hz or greater, more preferably1 Hz or greater). This enables the concept of driving a vehicle atnormal surface street speeds (e.g., 35 miles per hour) whileaccumulating useful gas concentration and wind measurement data.However, embodiments of the invention do not depend critically on thegas detection technology employed. Any gas concentration measurementtechnique capable of providing gas concentration measurements can beemployed in some embodiments.

Although the gas concentration measurements are preferably performedwhile the gas measurement device is moving, at least some gasconcentration measurements can be performed while the gas concentrationmeasurement device is stationary. Such stationary gas concentrationmeasurements may be useful for checking background gas concentrations,for example. While real-time measurements are preferred, post analysisof more sparsely sampled data, e.g., via vacuum flask sampling and lateranalysis via gas chromatography or other methods, may be used in someembodiments. Optionally, measurements can be made on different sides ofthe road or in different lanes to provide more precise localization ofthe leak source. Optionally, the present approaches can be used inconjunction with other conventional methods, such as visual inspectionand/or measurements with handheld meters to detect emitted constituents,to further refine the results. Optionally, measurements can be made atreduced speed, or with the vehicle parked near the source, to provideadditional information on location and/or source attribution.

Optionally, the system can include a source of atmosphericmeteorological information, especially wind direction, but also windspeed or atmospheric stability conditions, either on-board the vehicleor at a nearby location. The stability of the atmospheric conditions canbe estimated simply from the wind speed, the time of day, and the degreeof cloudiness, all of which are parameters that are available either onthe vehicle or from public weather databases. Optionally, the computersystem can include an on-board video camera and logging system that canbe used to reject potential sources on the basis of the local imagerycollected along with the gas concentration and wind data. For example, ameasured emissions spike could be discounted if a vehicle powered bynatural gas passed nearby during the measurements. Optionally, repeatedmeasurements of a single location can be made to provide furtherconfirmation (or rejection) of potential leaks. Accordingly, the scopeof the invention should be determined by the following claims and theirlegal equivalents.

What is claimed is:
 1. A method comprising: performing a gas leakdetection survey by measuring a plurality of gas concentration values ata corresponding plurality of measurement points using a gasconcentration measurement device carried by a moving vehicle along asurvey path, determining a plurality of geospatially-referencedlocations of the plurality of measurement points using a positioningdevice carried by the moving vehicle, and measuring a plurality of winddirection values using a wind direction determination device carried bythe moving vehicle, wherein the plurality of wind direction values aremeasured within a time interval during which the plurality of gasconcentration values are measured, the time interval being less than orequal to 2 minutes; determining, using at least one hardware processor,a range of search directions to a suspected gas leak source extendingfrom a location of a gas concentration measurement point according to avariability of the wind direction values, wherein the range of searchdirections defines a width of a graphical indicator of the variabilityof the wind direction values; and employing the range of searchdirections to identify a location of the suspected gas leak source, tofacilitate a repair of the gas leak source.
 2. The method of claim 1,wherein the variability comprises a variance or standard deviation ofthe wind direction values.
 3. The method of claim 2, wherein the rangeof search directions is proportional to the standard deviation of thewind direction values.
 4. The method of claim 1, further comprisingemploying the at least one processor to determine a representative winddirection relative to the gas concentration measurement point.
 5. Themethod of claim 4, wherein the representative wind direction comprises amean, median or mode of the wind direction values.
 6. The method ofclaim 1, further comprising employing the at least one processor todetermine, according to wind speed data relative to the gasconcentration measurement point, a maximum detection distance valuerepresentative of an estimated maximum distance from the suspected gasleak source at which a gas leak is detectable by the survey.
 7. Themethod of claim 6, wherein the maximum detection distance value isfurther determined according to data representative of atmosphericstability conditions along the survey path.
 8. The method of claim 1,wherein the location of the gas concentration measurement pointcorresponds to a peak in gas concentration identified from the pluralityof gas concentration values.
 9. An apparatus comprising: a gasconcentration measurement device configured to measure a plurality ofgas concentration values at a corresponding plurality of measurementpoints along a survey path travelled by a moving vehicle conducting agas leak detection survey, the vehicle carrying the gas concentrationmeasurement device; a positioning device configured to determine aplurality of geospatially-referenced locations of the plurality ofmeasurement points; a wind direction determination device configured tomeasure a plurality of wind direction values within a time intervalduring which the plurality of gas concentration values are measured, thetime interval being less than or equal to 2 minutes; and at least onehardware processor coupled to the gas concentration measurement device,the positioning device, and the wind direction measurement device, andprogrammed to determine a range of search directions to a suspected gasleak source extending from a location of a gas concentration measurementpoint according to a variability of the wind direction values, whereinthe range of search directions defines a width of a graphical indicatorof the variability of the wind direction values.
 10. The apparatus ofclaim 9, wherein the variability comprises a variance or standarddeviation of the wind direction values.
 11. The apparatus of claim 10,wherein the range of search directions is proportional to the standarddeviation of the wind direction values.
 12. The apparatus of claim 9,wherein the at least one processor is configured to determine arepresentative wind direction relative to the gas concentrationmeasurement point.
 13. The apparatus of claim 12, wherein therepresentative wind direction comprises a mean, median or mode of thewind direction values.
 14. The apparatus of claim 9, wherein the gasconcentration measurement point corresponds to a peak in gasconcentration identified from the plurality of gas concentrationmeasurements.
 15. The apparatus of claim 9, wherein the at least oneprocessor is further programmed to determine, according to wind speeddata relative to the gas concentration measurement point, a maximumdetection distance value representative of an estimated maximum distancefrom the suspected gas leak source at which a gas leak is detectable bythe survey.
 16. The apparatus of claim 15, wherein the maximum detectiondistance value is further determined according to data representative ofatmospheric stability conditions along the survey path.
 17. A methodcomprising: performing a gas leak detection survey by measuring aplurality of gas concentration values at a corresponding plurality ofmeasurement points using a gas concentration measurement device carriedby a moving vehicle along a survey path, determining a plurality ofgeospatially-referenced locations of the plurality of measurement pointsusing a positioning device carried by the moving vehicle, and measuringa plurality of wind direction values and a plurality of wind speedvalues using at least one wind measurement device carried by the movingvehicle, wherein the plurality of wind direction values and theplurality of wind speed values are measured within a time intervalduring which the plurality of gas concentration values are measured, thetime interval being less than or equal to 2 minutes; and determining,using at least one hardware processor, a probable direction and anestimated maximum distance to a suspected gas leak source relative to ageo-referenced location of a gas concentration measurement point,wherein the maximum distance represents an estimated maximum distancefrom the suspected gas leak source at which a gas leak is detectable bythe survey, wherein the probable direction is determined according tothe plurality of wind direction values, wherein the maximum distance isdetermined according to the plurality of wind speed values, and whereinthe maximum distance defines a length of a graphical indicator of themaximum distance.
 18. The method of claim 17, wherein the maximumdistance is further determined according to data representative ofatmospheric stability conditions in an area of the gas concentrationmeasurement point.
 19. The method of claim 17, wherein the location ofthe gas concentration measurement point corresponds to a peak in gasconcentration identified from the plurality of gas concentration values.20. An apparatus comprising: a gas concentration measurement deviceconfigured to measure a plurality of gas concentration values at acorresponding plurality of measurement points along a survey pathtravelled by a moving vehicle conducting a gas leak detection survey,the vehicle carrying the gas concentration measurement device; apositioning device configured to determine a plurality ofgeospatially-referenced locations of the plurality of measurementpoints; at least one wind measurement device configured to measure aplurality of wind direction values and a plurality of wind speed valueswithin a time interval during which the plurality of gas concentrationvalues are measured, the time interval being less than or equal to 2minutes; and at least one hardware processor coupled to the gasconcentration measurement device, the positioning device, and the atleast one wind measurement device, and programmed to determine aprobable direction and an estimated maximum distance to a suspected gasleak source relative to a geo-referenced location of a gas concentrationmeasurement point, wherein the maximum distance represents an estimatedmaximum distance from the suspected gas leak source at which a gas leakis detectable by the survey, wherein the probable direction isdetermined according to the plurality of wind direction values, whereinthe maximum distance is determined according to the plurality of windspeed values, and wherein the maximum di stance defines a length of agraphical indicator of the maximum distance.
 21. The apparatus of claim20, wherein the maximum distance is further determined according to datarepresentative of atmospheric stability conditions in an area of the gasconcentration measurement point.
 22. The apparatus of claim 20, whereinthe location of the gas concentration measurement point corresponds to apeak in gas concentration identified from the plurality of gasconcentration values.
 23. A non-transitory computer-readable mediumencoding instructions which, when executed by a computer system, causethe computer system to: receive a plurality of gas concentration valuesmeasured at a corresponding plurality of measurement points using a gasconcentration measurement device carried by a moving vehicle along asurvey path, receive a plurality of geospatially-referenced locations ofthe plurality of measurement points determined using a positioningdevice carried by the moving vehicle, and receive a plurality of winddirection values and a plurality of wind speed values measured using atleast one wind measurement device carried by the moving vehicle, whereinthe plurality of wind direction values and the plurality of wind speedvalues are measured within a time interval during which the plurality ofgas concentration values are measured, the time interval being less thanor equal to 2 minutes; and determine a probable direction and anestimated maximum distance to a suspected gas leak source relative to ageo-referenced location of a gas concentration measurement point,wherein the maximum distance represents an estimated maximum distancefrom the suspected gas leak source at which a gas leak is detectable bythe survey, wherein the probable direction is determined according tothe plurality of wind direction values, wherein the maximum distance isdetermined according to the plurality of wind speed values, and whereinthe distance defines a length of a graphical indicator of the maximumdistance.
 24. The medium of claim 23, wherein the maximum distance isfurther determined according to data representative of atmosphericstability conditions.
 25. A non-transitory computer-readable mediumencoding instructions which, when executed by a computer system, causethe computer system to: receive a plurality of gas concentration valuesmeasured at a corresponding plurality of measurement points using a gasconcentration measurement device carried by a moving vehicle along asurvey path; receive a plurality of geospatially-referenced locations ofthe plurality of measurement points determined using a positioningdevice carried by the moving vehicle, and receive a plurality of winddirection values measured using a wind direction determination devicecarried by the moving vehicle, wherein the plurality of wind directionvalues are measured within a time interval during which the plurality ofgas concentration values are measured, the time interval being less thanor equal to 2 minutes; and determine a range of search directions to asuspected gas leak source extending from a location of a gasconcentration measurement point according to a variability of the winddirection values, wherein the range of search directions defines a widthof a graphical indicator of the variability of the wind directionvalues.